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We compute the phase diagram of the HgTe/CdTe quantum wells in the 3 dimensional (3D) 
parameter space of Dirac mass, Fermi level and disorder strength. The phase diagram reveals the 
Quantum spin-Hall, the metallic and the normal insulating phases. The phase boundary of the 
Quantum spin-Hall state is shown to be strongly deformed by the disorder. Taking specific cuts 
into this 3D phase diagram, we recover the so called topological Anderson insulator (TAI) phase, 
but now we can demonstrate explicitly that TAI is not a distinct phase and instead it is part of the 
Quantum spin-Hall phase. The calculations are performed with -conserving and S'^^-nonconserving 
Hamiltonians. 

PACS numbers: 73.43.-f, 72.25.Hg, 73.61.Wp, 85.75.-d 



The HgTe/CdTe quantum wells were the first candi- 
dates as Quantum spin-Hall (QSH) insulators. Theory 
predicted that HgTe/CdTe quantum wells can be tuned 
to enter the QSH phase, where one should observe ro- 
bust conducting helical channels flowing around the edges 
of the samples.^ This prediction was partially confirmed 
by experiments,^ which showed quantized direct conduc- 
tance for small contact distances. The conductance mea- 
surements, however, consistently showed a decay of the 
direct conductance with the distance between the con- 
tacts, a hallmark of localization. This behavior spurred 
several studies that looked into the robustness of the 
edge modes of HgTe/CdTe quantum wells against de- 
phasing,^^ and strong disorder (references given below). 

While pursuing such program, one study detected a 
seemingly new phase, which was dubbed the topological 
Anderson insulating (TAI) phase. ^ Like the QSH phase, 
TAI was found to display perfectly conducting channels 
at the edges of the samples, in spite of the presence of 
strong disorder. The parameter region where this phase 
was observed seemed to be totally disconnected from the 
other phases, so one could legitimately conclude that TAI 
is indeed a distinct phase of the system. Specifically, 
TAI was observed at large disorder, for Dirac masses 
that would give a topologically trivial system if the dis- 
order were turned off, and for Fermi levels well inside 
the conduction band of the clean system. These findings 
were subsequently confirmed by independent studies,^ ^ 
thought sometime different interpretations were given. 

To be more concrete, let us introduce the relevant 
model. As shown in Ref.[TJ the "low energy" band theory 
of the clean HgTe/CdTe wells can be captured by a 2D 
effective Hamiltonian (written in the momentum space): 

/ h{k) r(fc) \ 
- r(fc)t h*{-k) j' 

where h{k)=e{k)-\-d{k)-a, with a={ax, Cy^ az) encoding 
the Pauli's matrices, and T{k) is a /S^-nonconserving in- 
teraction. The behavior of the measured energy bands 
in the proximity of the T-point can be captured by the 
following expression for d{k): 

d = (A sin kx^A sin ky^ A — 2B{2 — cos kx — cos ky)). (2) 



In the absence of disorder and with T{k) set to zero, Hq 
displays a topological phase if 0<A/5<4 and 4<A/5<8 
(the insulating gap closes when M/B=0^ 4 and 8), and 
a topologically trivial phase if A/B<0 or A/B>S. 

Working within the Born approximation, a second 
stud}/^ on disordered HgTe/CdTe quantum wells showed 
that disorder renormalizes the Dirac mass term A and 
the chemical potential of the system. The renormaliza- 
tion occurs in the "right" direction, so as to give rise 
to the TAI phase. Based on this analysis, this reference 
concluded that TAI occurs because of the real gap renor- 
malization. 

Numerically, the phase diagram of strongly disordered 
QSH insulators was computed in the {Ef^ W) plane {W 
= disorder strength) using the transfer matrix approach,!^ 
by probing the conductance of the edges, ^ ^ and by level 
statistics analysis and computation of the quantized bulk 
invariants.'^ The phase diagram was also computed in 
the (A, W) plane using the transfer matrix method.^' 
When computed in the {Ef^ W) plane, the TAI phase ap- 
pears to be disconnected (distinct) from the other phases, 
but when computed in the (A, W) plane, the TAI phase 
appears to be connected to, or be a part of the QSH 
phase. 

The studies mentioned above provided diverging views 
about the TAI phase, which led to a splitting of the com- 
munity following or working on this effect, with one side 
claiming that TAI is a novel phase and the other side say- 
ing that TAI is connected to the QSH phase, and there- 
fore is not a new phase. The issue can be settled down by 
a computation of the global phase diagram of the system 
in the 3D parameter space of {Ef^ A, W)^ in which case 
one can follow the evolution of the phase diagram in all 
three directions and therefore record the emergence (or 
the lack of it) of various disconnected phase components. 
This is exactly what we are set to do in this work. 

We have recently developed an efficient numerical 
method to evaluate the spin-Chern number in the pres- 
ence of disorder .I^^HI^ The methodology stems from a 
combination of the non-commutative theory of the spin- 
Chern number (applicable to iS^-nonconserving models 
too), ^ and a novel numerical algorithm to evaluate topo- 
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logical invariants in the real-space representationflSIill In 
the present work, we will use such computations to map 
the regions where the spin-Chern number displays quan- 
tized values and ultimately to construct the 3D phase di- 
agram of HgTe/CdTe QSH wells in t he prese nce of strong 
disorder. We have shown in the pastp^^^^^^ and we have 
verified this again, that such calculations are in perfect 
agreement with other methods of investigations, such as 
level statistics analysis. 

Based on our results, we can now show explicitly that 
the TAI phase is the same as the Quantum spin-Hall 
phase, and that the apparent distinct character of the 
TAI phase seen in the original papei^ is due to the fact 
that the phase diagram was examined only along a 2D 
section of this 3D phase diagram. We will show exactly 
where that section occurs in our 3D phase diagram. In 
fact, we will show that our 3D phase diagrams are in good 
quantitative agreement with all the previous numerical 
calculations. 

Besides resolving the nature of the TAI phase, the 
present work intends to showcase the efficiency and ac- 
curacy of the calculations based on the non- commutative 
spin-Chern number. 

I. THE MODEL SYSTEMS 

The computations will be carried out for models de- 
scribed by Hq plus disorder. The generic form of the 
diagonal term is 

e{k) = C — 2D{2 — cos kx — cos ky)^ (3) 

and the minimal form of the matrix T{k) was derived in 
Ref.El 

T(k) = iA( sin/cx-^sin/c^ \ 

^ ^ y sin kx ^ i sin ky J ' ^ ^ 

The real space representation of Hq can be constructed 
on a square lattice where each vertex n carries four quan- 
tum states |n,a,cr), where a = ±1 (= isospin) labels 
the s or the p character of the bands and a = ±1 the 
spin up and down configurations. On the Hilbert space 
spanned by |n,a,cr), the translational invariant Hamil- 
tonian takes the form (a = lattice spacing): 

Ho = §{ti,o + + ^0,1 + ^0,-1 - 4 + C/D) 

(5) 

+ ^0^(^1,0 + ^-1,0 + ^0,1 + ^0,-1 - 4 + M/B) 
where t 

m.ki ^ 1 ^1 s^nd T(j are the translations, spin, 
isospin and fiipping operators defined below: 

tm,k\ni, 712, a, cf) = \ni + m, n2 + k,a,a), 
(j|n, a, <j) = cr|n, a, <j), a|n, a, a) = a|n, a, <j), (6) 
rcr|n,Qf,cr) = |n,a, — <j), ra\Ti^a^(j) = |n, — a,cr). 



When the parameters are given the values f^^ A = 364.5 
meV nm, B = -686 meV nm^ C = 0, D = -512 meV 
nm^, A = and a = 5 nm, the Hamiltonian Hq accu- 
rately reproduces the relevant band structure of the prop- 
erly tuned HgTe/CdTe wells in the QSH regime. Our 
calculations will be carried with these realistic parame- 
ter values but also with the theoretical values: A = 1, 
5 = 1, C = 0, D = 0, A = (and A 0.5), and a set 
to one. These later theoretical values have been used in 
Ref. i8| which contains one of the most accurate calcula- 
tions to date for disordered topological insulators. For 
this reason, we decided to use Ref. [SI for comparison and 
for assessing the accuracy of our calculations. 

The computations with disorder will be done with the 
Hamiltonian: 

= Ho^W ^ ujn,a\n,a,a){n,a,a\, (7) 

n,a,a 

where ujn (X are random amplitudes uniformly distributed 
in the interval [— ^, This is just a crude approxima- 
tion of the disorder in QSH wells, where the disorder is 
probably mostly due to the random displacements of the 
atoms from the perfect crystalline structure. The lead- 
ing physical effect of such displacements will be a change 
of the overlap integrals of the atomic orbit als. Conse- 
quently, the disorder in QSH wells will primarily occur 
in the hopping amplitudes. We are in the process of in- 
vestigating such effects, but for now, we will follow the 
previous studies and work with the disordered Hamilto- 
nian of Eq. [7) 



II. THE NON-COMMUTATIVE SPIN-CHERN 
INVARIANT 

This invariant was discussed extensively in Refs. [TqUT21 
and here we will give only a brief account of it. Let 
be the operator a^ln, a, cr)=(j|n, a, a) and suppose that 
an exact diagonalization was performed for i^^. Given 
a Fermi level, one can compute the projector P^^ onto 
the energy spectrum below Ep. Furthermore, one can 
diagonalize the operator P^cFzPuj and what he will find is 
an eigenvalue spectrum that is symmetric relative to the 
origin and contained in the interval [—1,1]. We denote 
by P^ the spectral projector onto the positive/negative 
spectrum of P^azPu- Effectively, these projectors split 
the space of the occupied states into spin up and spin 
down sectors, and no 5';^ conservation is required by the 
procedure. Now, the projectors P^ fit into the non- 
commutative theory of the Chern invariant,!^ so one can 
define the non-commutative Chern numbers: 

C± = 2^z(tro{P^ [ - P±], 2, P^]] }>, (8) 

where the outer angular parentheses signify disorder av- 
erage, and tro is the trace over the states at site n = 0. 
Also, X is the position operator. The non-commutative 
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FIG. 1. (Color online) The numerical values of the spin-Chern 
number Cs as function of A, when the other parameters were 
fixed at: A = B = a = 1, A = 0, = 5.71 and Ef = 0. 
Cs was computed for three lattices of increasing size: 30x30 
(triangles), 40x40 (squares) and 50x50 (circles). The shaded 
region indicates the region where \Cs\ > 0.5, which is being 
identified with the QSH phase. 



spin-Chern number is defined as 



C, = -(C+-C_). 



(9) 



An immediate consequence from Ref. 16 is that Cs 
takes quantized values as long as the spectrum of P^&zPoj 
is localized near the origin so that the matrix elements 
of decay sufficiently fast, more precisely, as long as: 



ma,a\P^ 



< oo. 



(10) 



One can rigorously show that the region where IC^I = 1 is 
inside the QSH phase, and that the region where = 
is inside the normal insulating phaseP^ Further investi- 
gations have shown that the delocalization of P^ occurs 
simultaneously with the delocalization of the full projec- 
tor P^\^^^ The reason for this is because the mobility 
gap of P^&zPoj is insensitive to spin-independent disor- 
der, so its eigenstates near the origin remain localized as 
long as P^ stays localized. The practical consequence of 
all this is that the region where Cg takes the quantized 
values ±1 actually coincides with the QSH phase region 
and the region where Cs takes the value coincides with 
the normal insulating phase. The region where Cs takes 
non-quantized values can be identified with the metallic 
phase. These conclusions have been verified by extensive 
numerical computations.^^ 

The non-commutative formulas of Eq. [S] makes sense 
only in the thermodynamic limit but in practice we can 
only work with finite size systems. For a finite N x N 
square lattice with periodic boundary conditions, the fol- 
lowing formula was derived 
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TV{P± [ - i [xi , P±J ,-ilx2,P^\]}, (11) 
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FIG. 2. Global view of all our numerical data for A 





(see Section III A), obtained on the 50x50 lattice. A, B and 
a took the same values as in Fig. [l] Cs is represented by a 
surface plot, with Ef and A on the horizontal axes and Cs on 
the vertical axis. There is one surface plot for each W value 
(shown on the right) we considered in our study. The surface 
plots are shown from three angles: sideway from Ef direction 
(first column), sideway from A direction (second column) and 
from above (third column). 



where the trace is over all the states, and 



(12) 



Above, Q is taken of the order of N/2 and the c^'s are 
solutions of the following linear system of equations: 
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With these choices, the finite size formulas from Eq. [TT] 
converges exponentially fast to the exact formulas given 
in Eq. |8j as long as the Fermi level is in a mobility gap. 
The exponential convergence slows down when the Fermi 
level nears a mobility edge. Consequently, the most chal- 
lenging part of the calculations is to achieve convergence 
near the mobility edges, same as to say, near the phase 
separation lines. 



From a numerical point of view, Eq. 11 is straightfor- 
ward to implement. Indeed, note that the matrix ele- 
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FIG. 3. The contour levels Cs — 0.5 obtained from the surface 
plots of Fig. [2] and from similar data obtained on the 30x30 
and 40x40 lattices. When the calculations are converged with 
the size of the lattice, the region inside these contours contains 
the QSH phase. 



merits of [^i,^^^] take the simple form: 



(n,a,cr|[fi,P^J|n',a',a') = £ c^x 

m=l (14) 



sm 



Also note that are self- averaging, so the uj- 

dependence becomes weaker and weaker as the volume is 
taken larger and larger. In practice, we still need to aver- 
age over a few disorder configurations in order to obtain 
relatively smooth data. To summarize, the computation 
will proceed as follows: 

• Chose a value for A, Ep and W . 

• Generate a random potential and diagonalize the 
resulting R^^. 

• Construct the projector P^^ onto the states with 
energy below Ep. 

• Construct and diagonalize the operator P^azPuj- 

• Construct the projectors P^ onto the posi- 
tive/negative spectrum of P^cTzPuj- 



• Evaluate the finite size formulas Eq. 11 



Compute the spin-Chern number using Eq. |9] 
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FIG. 4. The 3D phase diagram of the model with A = 0, as 
derived from the contours of Fig.|3] The 3D region delimitated 
by the contours and the dotted line contains the QSH phase. 
The shaded region indicates the region which contains the 
so called topological Anderson insulating phase. The dotted 
lines were obtain by extrapolation. 



• Analyze the spin-Chern number data. 

We end this section by mentioning that, for time- 
reversal symmetric models, the spin-Chern number mod- 
ulo 2 coincides^^ with the Z2 invariant introduced by 
Kane, Mele and Fu in Refs. ITtI and ITSl For example, the 
analytic calculations of Re f-TO" showed explicitly that, for 
a certain ^'^-nonconserving model in the clean limit, the 
kernel appearing in the computation of the spin-Chern 
number is identical to the Pfaffian function used in the 
computations of the Z2 invariant (in the spirit of Ref. [T7|) 
for the same model, performed in Refs. 20 and 2T, For 
this analytically solvable model, the Z2 and the spin- 
Chern invariants were explicitly shown to give identical 
phase diagrams. 



III. MODEL CALCULATIONS AND ANALYSIS 

In this section we set the parameters to the theoretical 
values mentioned above: A = B = a = 1 and zero in 
rest, except for A which will be given the values and 
0.5. When doing so, the model becomes particle- hole 
symmetric and the phase diagram becomes mirror sym- 
metric relative to the plane A = 4. This is of great help 
when computing the 3D phase diagram, since we only 
need to cover the A < 4 and Ep < region (one quar- 
ter of the whole parameter space). This is what we will 
do and, once the calculations are completed, we recon- 
struct the full phase diagram by using these mentioned 
symmetries. 
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FIG. 5. (Color online) The numerical values of the spin-Chern 
number Cs as function of A, when the other parameters were 
fixed at: A = 5 = a = 1, A = 0.5, W = 5.71 and Ef = 0. 
Cs was computed for three lattices of increasing size: 30x30 
(triangles), 40x40 (squares) and 50x50 (circles). The various 
shaded regions indicate the regions where a) 1.5 < \Cs\ < 
0.9, which is being identified with the metallic phase and b) 
\Cs\ > 0.9, which is being identified with the QSH phase. 



A. -conserving calculations 

We set here A to zero, in which case Sz commutes with 
the Hamiltonian and the model decouples into two Chern 
insulators, one on the spin- up sector and the other one 
on the spin-down sector. The model now belongs to the 
unitary class so one expects to see a sharp transition be- 
tween the QSH and normal insulating phases-f^HIS^ The 
states strictly between the QSH and normal insulating 
phase are necessarily delocalized.^^ 

We first look in more detail at the numerical results. 
For this, we fix W — 5.71 and Ep = ^ and examine the 
values of Cs as A is varied from —2 to 10. Graphs of 
these values are shown in Fig. [l] for lattices of increasing 
sizes: 30x30, 40x40 and 50x50. The level of disorder 
is high, the value W = 5.71 being comparable with the 
width of the clean bands and much larger than the occur- 
ring insulating gaps. Consequently, the insulating gaps 
are completely filled with dense insulating spectrum in 
all calculations presented in Fig. [l] Cs was averaged 
over four disordered configurations. We have not consid- 
ered larger sizes or additional disordered configurations 
because the resulting phase diagram is already well con- 
verged and the averaged Cs appears smooth. As one can 
see from Fig. [ij there are regions where Cs takes quan- 
tized values of and ±1, but also regions where Cs takes 
non-quantized values. The regions with Cs = can be 
safely regarded as corresponding to the normal insulat- 
ing phase, while those with = 1 as corresponding to 
the QSH phase. It is known analyt ically that |Cs|=l im- 
plies existence of edge states.'^^^ To resolve the region 
where Cs takes non-quantized values, which occurs pre- 
cisely near the phase boundaries as discussed above, we 
relay on the behavior of Cs when increasing the size of 
the lattice. By examining each point of the diagram, we 
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FIG. 6. Global view of all our numerical data for A = 0.5 



(see Section IIIB), obtained on the 50x50 lattice. A, B and 
a took the same values as in Fig. |5] Cs is represented by a 
surface plot, with Ef and A on the horizontal axes. There 
is one surface plot for each W value (shown on the right) we 
considered in our study. The surface plots are shown from 
three angles: sideway from Ef direction (first column), side- 
way from A direction (second column) and from above (third 
column) . 



see that, for specific values of A, Cs moves towards ±1, 
and, for other values of A, Cs moves towards 0. Thus, 
the size dependence of Cs gives us a practical method to 
resolve the regions where the convergence is slow. 

As an empirical fact, we found that the points 
above/below the values ±0.5 move up/down when the 
size of the lattice is increased, while the other points move 
down/up. Therefore, with a good measure, we can take 
the points where |Cs |=0.5 as the phase boundary between 
the QSH and normal insulating phases. This criterion is 
exemplified on the data shown in Fig.[l] which then leads 
to the ID phase diagram shown in the same figure. 

Fig. [2] gives a global view of the data obtained on the 
50x50 lattice, by plotting Cs as a function of both A and 
Ef for several disorder strengths W. All our conclusions 
in this section are based on this data. From Fig. [2j the 
reader can get a sense of how sharp are the regions where 
Cs takes quantized values, how precise this quantization 
is (it is very precise), and what is the extent of the regions 
where the calculations are not fully converged. From this 
global mapping we can extract the 3D phase diagram of 
the model. To do that, for each surface plot in Fig. [2j 
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we draw the contours corresponding to the level values 
|Cs 1=0.5 and the results are shown in Fig. [s] We have 
placed the results for the different lattice sizes near each 
other so that one can examine the convergence of the 
phase diagram. It is quite evident that the contours are 
well converged with the size of the system for W values 
up to 7.14. The computed QSH phase still has a small 
growth for W=S.57 as the size is increased. But overall, 
we believe the 50x50 lattice calculations gives a fairly 
well converged 3D phase diagram of the model, which is 
shown in Fig. [4j 

The phase diagram shown in Fig. [4] is in excellent 
quantitative agreement with the data reported in Ref. |8] 
(where only the section Ep = was examined). As al- 
ready noted in Ref. 8, the phase boundary is strongly 
reshaped by disorder; the QSH phase is monotonically 
downsizing and drifts away from the A = 4 plane as W 
is increased. Because of these particularities, the QSH 
phase extends in the regions A < and A > 8, which 
correspond to the normal insulating phase when 1^=0. 
We dubbed this region the TAI region, because here is 
where the TAI phase will be observed if one would ana- 
lyze only a slice (for example A=-0.5) of the 3D phase 
diagram. But now it becames clear that the TAI phase 
is not a new phase but is part of the QSH phase, whose 
phase boundary was strongly reshaped by the disorder. 
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B. AS^-nonconserving calculations 

Here we set the parameters at A=B=a=l and A=0.5, 
so that a large Sz non- conserving potential is present. In 
this case, the model belongs to the symplectic class, so 
one expects to see a metallic phase in between the QSH 
and normal insulating phases. The task of computing 
the phase diagram is much more difficult now, but never- 
theless can be accomplished using the spin-Chern num- 
ber. 

Let us explore the numerical data in detail. As before, 
we fix two parameters: 1^=5.71 and £^^=0, and then 
examine the values of Cs as A is varied from —2 to 10. 
The values are graphed in Fig. [5j for three lattice sizes: 
30x30, 40x40 and 50x50. Again, one can see regions 
where Cs takes quantized values ±1 and 0, which can 
be safely regarded as belonging to the QSH and normal 
insulating phases, respectively. There are also regions 
where Cs doesn't take quantized values. The metallic 
phase is contained within this region, but of course there 
will be an inherent uncertainty in establishing the exact 
boundaries of the phases. 

Our criterion for phase delimitation will be as fol- 
lows. We identify the QSH phase with the region where 
\Cs\ > 0.9 and the normal insulating phase with the re- 
gion where |Cs| < 0.15. The region where 0.15 < \Cs\ < 
0.9 will be identified with metallic phase. If we apply 
this criterion to the data in Fig. [5) we obtain the ID 
phase diagram shown on top of the graphs in the same 
figure. The value |Cs|=0.15, chosen to define the bound- 



FIG. 7. The contour levels \Cs\ = 0.15 and = 0.9 ob- 
tained from the surface plots of Fig. [6] and from similar data 
obtained on the 30x30 and 40x40 lattices. The contours 
are easy to distinguish because the 0.9 level contours are al- 
ways inside the 0.15 level contours. The regions inside the 
0.9 level contours give an approximate representation of the 
QSH phase, and the regions between the 0.15 and 0.9 level 
contours give an approximate representation of the metallic 
phase. 



ary between the metallic and normal insulating phase, 
may appear high but this was the lowest value we could 
consider and still obtain smooth boundary separations. 

In Fig. [6] we give a global view of the data obtained 
with a 50x50 lattice. All our conclusions for this sec- 
tion are based on this data. It is quite obvious that the 
transition from the QSH phase to the normal insulating 
phase is much blurrier this time (compare with Fig. [2]) , 
owing to the emergence of the metallic phase when A 7^ 0. 
We now apply the above criterion and draw the contours 
corresponding to |Cs|=0.15 and 0.9 for each independent 
surface plot shown in Fig. [g] (and similar plots generated 
with 30x30 and 40x40 lattices). The results are shown 
m Fig.[7l Here we again placed near each other the re- 
sults for the three different lattice sizes: 30x30, 40x40 
and 50x50, so that one can judge the convergence of the 
phase boundaries with the size of the system. In our opin- 
ion, the contours in Fig. [T] are well converged so we can 
proceed with the drawing of the phase diagram, which is 
shown in Fig. [Sj Our phase diagram is again in excellent 
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FIG. 8. The 3D phase diagram of the model with A = 0.5, as 
derived from the contours of Fig. [71 The 3D region dehmitated 
by the inside contours and the dotted hne contains the QSH 
phase. The 3D region between the contours and the dotted 
Unes Unes contains the metaUic phase. The shaded region 
indicates the region which contains the so called topological 
Anderson insulating phase. The dotted lines were obtained 
by extrapolation. 



1D0 
d 

100 




30x30 



40x40 



50x50 



Ef 






I 








Ef 








r 


A 








A 








A 



t 

274.4 



'100 -50 100 -SO 100 -so 



I 


:f 




Ef 




Ef 










A 




A 








^A 



200 
100 





?00 
100 


200 
TOO 




Ef 




100 -50 100 -50 100 -50 ^ 



Ef 


Ef 


Ef 









-100 -50 D 100 -50 m -50 




Ef 



Ef 



Ef 



Ef 



117.6 



78.40 



■100 -50 100 -50 100 -50 



-100 -SO 100 -SO 100 -SO 



Ef 


Ef 


Ef 









A I 



■100 '50 100 '50 100 -SO 



quantitative agreement with the one reported in Ref. El 
Examining this 3D phase diagram, we observe the 
same trends seen in the previous calculations, with the 
QSH phase monotonically downsizing and moving away 
from the plane A = 4. The QSH phase region is smaller 
now when compared with the case A = 0, but still there 
is a TAI region. In fact, as already mentioned in Ref. IH 
there is a region of A values where, by just increasing 
the disorder in the system, one will observe a sequence 
of phase changes from normal insulator to metal, then to 
QSH insulator, then to metal and back to normal insu- 
lator. 



IV. QUANTUM WELLS IN THE SPIN-HALL 
REGIME 

Here we set the parameters so that the band structure 
of Ho matches that of the HgTe/CdTe quantum well in 
the spin-Hall regime: 74=364.5 meV nm, 5=-686 meV 
nm^, (7=0, D=-512 meV nm^, A=0 and a=5 nm. These 
same parameters were used in the theoretical studies of 
Refs.'5'and^, which we will compare with. Note that the 
5'2-nonconserving potential is being turned off in these 
calculations. 

It was shown in these two mentioned references that 
the QSH phase region (and the one referred to as TAI) 
moves up in energy with the increase of the disorder am- 
plitude. No such behavior was observed in our previous 



FIG. 9. The contour levels \Cs\ = 0.9 obtained for the 
HgTe/CdTe QSH wlls computed on 30x30 and 40x40 lattices 
and for various disorder strenghts. When the calculations are 
converged with the size of the lattice, the region inside these 
contours contains the QSH phase. 

calculations, so we must conclude that this behavior is 
triggered by the presence of the diagonal term e{k). The 
self-consistent Born approximation analysis reported in 
Ref. 6 does not capture this fact; for example, Eq. 5 from 
this reference, giving the renormalized mass and chemi- 
cal potential, implies that the lower edge of the mobility 
gap moves down when increasing W. Another aspect re- 
vealed by Refs. 5 and 6 is that the TAI region moves 
deeper and deeper into the normal insulating region of 
the clean limit when W is increased. We will pay partic- 
ular attention to these two aspects. 

We have repeated the calculation and analysis of the 
spin-Chern number along the steps detailed in the pre- 
vious sections. In Fig. [9) we show the phase separation 
lines between the normal and the QSH phase. The crite- 
rion adopted here is that QSH phase corresponds to the 
regions where IC^I > 0.9. Note that this is different from 
the criterion used in our model calculation, and there is 
a good reason for this. Because the QSH phase is be- 
ing pushed at higher energies, the localization lengths of 
the localized states are much larger in the present cal- 
culations than those observed in the model calculations. 
As a result, the transition between the QSH and normal 
insulating phase is not as sharp as in the model calcula- 
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FIG. 10. The solid dots and solid lines represent 2D slices of 
the 3D phase diagram shown in Fig. [s] (corresponding to the 
50x50 lattice), taken at: (a) A = -13.5 meV and (b) A = 5.75 
meV. The dotted contour in (a) represents the boundary of 
the QSH phase, which we traced out from Fig. 2c of Ref. 
The dotted contour in (b) represents the boundary of the TAI 
phase, which we traced out from Fig. 2f of Ref. [5] 



tions. By choosing the high value of 0.9, we tried to make 
sure that the computed QSH phase is actually inside the 
exact QSH phase. In this way, we can be sure that, when 
the approximately computed QSH phase spills out into 
the normal phase, the exact QSH phase does the same 
thing. 

Fig. [9] reveals that no disconnected phase regions are 
emerging when increasing W. Instead, the boundary be- 
tween the QSH and normal phases moves continuously 
deeper and deeper inside the A > region. In the same 
time, the boundary moves up in energy. This behavior 
is very similar to the behavior of TAI phase mentioned 
above. Refs.[5]and|6]drew phase diagrams in the (VF, Ep) 
plane corresponding to A = —10 meV and A = 1 meV. 
These diagrams can now be viewed as slices of the 3D 
phase diagram shown in Fig. [9) taken at the appropriate 
planes A = — 10 and 1 meV. We want to show explicitly 
that this is the case. Because of the way we discretized A 
in our numerical calculations, we cannot consider these 
same A values, and instead we will draw the sections 
corresponding to A = —13.5 meV and A = 5.75 meV. 
These sections are shown in Fig. [Toj together with the 
QSH and TAI phase boundaries from Refs.^and^^ The 



quantitative agreement seen in Fig. [10] leaves no doubt 
that the phase diagrams reported in Refs. 5 and 6 are 
just particular slices of the 3D phase diagram reported 
in Fig. [9| 

Nevertheless, the results show that a material can be- 
come a QSH insulator by just increasing the disorder of 
the crystalline structure. It is a well established fact that 
a trivial Anderson insulator will eventually emerge at 
very large disorder strengths^^ so, at least in this regime, 
one can expect a monotonic shrinkage of the QSH part 
of the phase diagram as the disorder is increased. But 



at lower disorder strengths, we now know that the phase 
diagram can display complex trends. For the particular 
model of the HgTe wells and the particular disorder con- 
sidered in our study, we actually see an initial expansion 
of the QSH phase, accompanied by a drift. For weak dis- 
order strengths, one can use a perturbative calculation to 
analyze the trends, and now it is well known from Ref. 6] 
that the first correction to the Dirac mass has a positive 
sign, hence a drift to right of the QSH phase is expected. 
What was also interesting in Ref. 6 was that the pertur- 
bative analysis remained accurate up to quite large dis- 
order strengths, so in some sense the general trends are 
set by the behavior at very weak disorder. It remains, 
however, to be explored how universal these conclusions 
are. For example it will be interesting to see what hap- 
pens when a different kind of disorder is considered, like 
the one affecting the hopping terms as discussed in one 
of our sections. 



V. CONCLUSIONS 

We have computed the phase diagram of a Quantum 
spin-Hall model Hamiltonian in the 3D parameter space 
of Dirac mass, Fermi energy and disorder strength. The 
analysis was based on the computation of the spin-Chern 
number, which was shown to display quantized values in 
the presence of large disorder and large ^^-nonconserving 
interactions. Working with certain theoretical values, 
we showed first that our computations are in excellent 
quantitative agreement with the 2D phase diagrams re- 
ported in Ref. '8', The 3D phase diagrams confirm the ab- 
sence of disconnected phase pockets and instead reveal a 
strong disorder-induced deformation of the phase bound- 
aries. Secondly, working with parameters appropriate for 
a HgTe/CdTe quantum well tuned for the quantum spin- 
Hall effect, we observed a similar disorder-induced de- 
formation of the phase boundaries and no disconnected 
phase pockets. We have shown that the 2D phase di- 
agrams reported in Refs. 5 and 6 correspond to slices 
appropriately taken from our 3D phase diagram. This 
shows explicitly that the so called TAI phase is not a 
new and distinct phase. Instead, TAI is part of the QSH 
phase whose boundary is strongly deformed as the disor- 
der amplitude is increased. 
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